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SPECTRAL  LINE  INTERFEROMETRY 
AND  INTERFEROMETER  NOISE  ANALYSIS 


I.  INTRODUCTION 

1 

The  detection  of  18 -cm  radiation  from  interstellar  OH  molecules  has  shown  the  need  for 

adding  the  frequency  dimension  to  interferometry.  Because  of  the  inability  of  large  parabolic 

antennas,  such  as  the  NHAO  14  0-foot  antenna  and  the  210- foot  antenna  in  Australia,  to  resolve 

the  18-cm  spectral  lines,  several  spectral  line  interferometers  have  been  built  to  measure  the 

2-4 

spatial  distribution  of  OH  emission.  Basically,  spectral  line  interferometry  is  no  different 
from  conventional  interferometry  in  which  no  frequency  analysis  is  performed;  however,  the 
spectral  line  interferometer  is  more  versatile  than  the  more  conventional  interferometer  and 
lias  led  to  some  improvements  in  conventional  interferometry.  In  particular,  the  high-speed 
computers  and  digital  processors  used  in  spectral  analysis  are  easily  adapted  to  perform  fringe 
rotation  and  delay  compensation  as  well  as  the  spectral  analysis. 


II.  PRINCIPLES  OF  INTERFEROMETRY 

Fourier  theory  relates  the  time  domain  electric  field  E(r,  t)  to  its  Fourier  components 

-*  -*■ 

E(r,  co)  in  the  frequency  domain  by  the  Fourier  integrals 


r:,  -  f  ,'V*  v  icot  dco 

Mr,  t)  ^  h(r,  w)  e 

E(r,  a.')  ^  K(r,  t)  e  dt 


(11-1) 


(11-2) 


In  an  exactly  analogous  wav,  the  Fourier  component  E(r,  co)  can  be  decomposed  into  spatial  com¬ 
ponents  E(k,  co)  where  each  spatial  component  is  a  plane  wave 


so  that 


and 


He  (E(k,  co )  e 


icot-ik-r 


} 


r> n  .,-*■*  dk  dk 

E(?.w)  =  V  D(k,c)e-lk-r 

(2tt) 


E(k,  to)  =  j  J  E(r,  to)  elk’r  drxdi 


(11-3) 


(11-4) 


(II-5) 
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Integration  over  only  two  dimensions  is  necessary  as 

|k|  =  ^  (II  -6) 

thus,  only  two  components  of  k  are  independent. 

Theorems  of  the  Fourier  theory  of  time  functions  carry  over  into  the  two-dimensional  Fourier 
theory  of  plane  waves.  The  relation  of  the  power  spectral  density  S(r,  o;)  to  the  autocorrelation 
function  <E(r,  t)  x  II* ( r,  t  —  t )^> 


S(r,  w)  -  E(r,  w)  x  H*(r,  w)  -  \  <E(r,  t)  x  H*(r,  t  -  t)>  e 


I 


dr 


(H-7) 


which  holds  for  a  time  function  assumed  zero  outside  the  interval 
O^t^T 

has  a  two-dimensional  equivalent 


S(k,  w)  =  E(ic,  w)x  H*(k,  u)  =  \\  <E(r,  u)  x  H*(r  -  y,  «)>  ,  elk‘Y  dy  dy  (11-8) 

JJ  H  x  y 

and  inverse  transform 

-*  pe  ^  .7*  dk  dk 

<E(r,  w)  x  Il*(r  -  y,  w)>  \  \  S(k,  w)  e‘lK'Y  —  -J-  (11-9) 

R  JJ  (2-) 


where  <  >  ^  denotes  a  time  average 


E(r,  t)  X  H(r,  t  -  t) 
T 


dt 


(11-10) 


and  ^  X*  a  spatial  average.  Thus,  the  power  in  the  direction  k  is  given  by  the  two-dimensional 
R 

transform  of  the  spatial  crosscorrelation  function.  The'  power  in  the  direction  k  per  unit  solid 
angle  and  frequency  interval  is  more  commonly  known  in  radio  astronomy  as  the  brightness 
distribution,  while  the  spatial  crosscorrelation  is  often  known  as  the  complex  visibility  function 
or  "fringe”  amplitude  and  phase.  When  two  antennas  are  used  to  measure  the  spatial  cross- 
correlation,  the  electric  fields  have  been  filtered  by  the  antenna  apertures.  If  the  two  antennas 
decompose  the  electric  field  vector  into  two  orthogonal  polarizations  a  and  b,  then 


T 


a,b 


(kx' 


k  ,w) 

y 


If  Xa^)Y^)elk'y  dyxdyy 

(g  (k  ,  k  ,  ic)  G*  ( k  ,  k  ,  u))\ 
xa  X  y  yb  X  y 


(11-11) 
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where  X(w)  and  Y  ( a; )  are  the  Fourier  components  of  the  antenna  outputs.  G  J  and  |  G  J  are 

x  y 

the  antenna  power  gains,  and  y  is  the  vector  from  the  phase  center  of  antenna  Y  to  the  phase 

center  of  antenna  X.  T  ,  is  the  brightness  temperature  for  the  Stokes  parameter  a,  b.  k  •  y/u; 

a ,  D 

is  the  delay  which  a  signal  wave  vector  k  will  undergo  in  reaching  the  phase  center  of  antenna 

X.  |k|  Yx/2?r  and  |k|  y y/27r  are  the  projected  baseline  components  in  wavelengths  which  produce 

fringe  spacing  of  Z 7r/  |  k |  y  r  and  2 7r/ )  k |  y  r  radians  in  the  x  and  y  directions.  The  sampling 

x  y 
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theorem  which  is  usually  applied  to  band-limited  functions  can  be  applied  to  the  complex  visibility 

and  its  transform.  When  the  source  distribution  is  limited  to  ±0  and  ±0  ,  then  the  complex 

x  _a  y  ~ \  5 

visibility  function  need  only  be  sampled  at  integral  multiples  of  (20  )_1  and  (20  )  wavelengths. 

x  y 

Also,  since  the  brightness  temperature  is  real 


(X  (w)  Y*(w)]  =  [X  (u)  Y*(w)]* 

1  a'  bv  'yx,Yy  1  av  bv  '-yx.-y 


(11-12) 


so  that  only  half  the  complex  visibility  plane  has  to  be  mapped. 

In  practice,  it  is  not  easy  to  sample  the  complex  visibility  plane,  as  it  is  physically  difficult 
to  move  antennas;  consequently,  in  many  interferometers  the  only  points  on  the  complex  visibility 
plane  that  can  be  sampled  are  those  on  the  contour  mapped  out  by  the  source  motion.  Sources 
with  fixed  right  ascension  and  declination  trace  out  an  ellipse  on  the  complex  visibility  plane  — 
a  portion  of  which  will  be  inaccessible  if  the  source  is  not  circumpolar.  If  azimuth  and  elevation 

A  A  A 

are  used,  coordinates  refer  to  the  Zenith  iy,  North  i  and  East  i^,  then 


_►  A  A  A 

y  =  D  sinE^i^  +  D  cosEn  sinA^i^  +  D  cosEn  cosA^L, 
H  Zi  H  is  L  is  is  IN 

-*  Cjl'  r  A  A  A 

k  =  —  {sinE^i^  +  cos  Ec  sinA0i^  +  cosEc  cosAeL-} 
c  1  SZ  S  SE  S  SN 


(11-13) 

(11-14) 


where  D  is  the  baseline  length,  E^  and  AR  are  the  elevation  and  azimuth  of  the  baseline,  and 
Ag  and  Eg  are  the  apparent  azimuth  and  elevation  of  the  source.  C  is  the  velocity  of  light  along 
the  baseline.  The  axes  of  rotation  on  some  antenna  mounts  do  not  intersect,  so  the  baseline 
becomes  dependent  on  the  source  position.  In  these  cases,  additional  vectors  have  to  be  intro¬ 
duced.  From  Eqs.  (11-13)  and  (11-14), 


intf 


(sinEB  sinEg  +  cos  E^  cos  Eg  cos(Ag 


■  ab)} 


(11-15) 


It  is  often  more  convenient  to  use  coordinates  centered  on  a  reference  position  in  the  sky  close 
to  the  source  being  studied,  in  which  case 


y  •  k 


intf 


{sin  6  sinS^  +  cos  6  cos  6  cos  (  Lc  —  LD)}  ( II  - 1 6 ) 

C  Ho  Ho  o  H 


where  L  is  the  hour  angle,  and  <5  is  the  declination.  In  this  coordinate  system,  k^  and  k^  can 

be  expressed  to  first  order  in  terms  of  the  angular  displacements  in  0^  right  ascension  and  0^ 

declination.  The  components  of  projected  baseline  in  the  direction  of  increasing  right  ascension 

and  declination  are  y  and  y  ,  where 

x  T  y 


yx=  D  cos6b  sin(Ls-LB) 

y  =  D{sin<5  cos  6  —  cos  <5  sin<5c  cos  (  Lc  —  L  )} 

y  Ho  Ho  o  H 


(H-17) 

(11-18) 


k  = 

x 


OJ 

C 


0 

X 


(11-19) 


k 


y 


C£ 

c 


0 

y 


(11-20) 


III.  METHODS  OF  CROSS  CORRELATION 

Although  the  spectral  visibility  function  Xa(w)  Y*(u;)  can,  in  principle,  be  determined  by 
Fourier  analysis  and  cross  multiplication  at  the  signal  frequency,  in  practice  the  signal  band 
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is  reduced  to  video  by  single-  or  multiple-frequency  conversion  stages.  After  mixing  with  a 
local  oscillator  of  frequency  taQ  and  phase  9,  the  Fourier  components  of  the  output  are 


x(w)  =  X(w  +  CO  )  e'10  +  X(a  -  co  )  e10 


(III-l) 


Thus,  for  an  upper-sideband  system 
X(ca)  =  x(u’  —  a'^)  e^ 


ca  >  0 


( III-  2) 


and 


iA^  t  . 

Y(ta)  =  y(ca  -  ic^)  e  c 


a*  >  0 


(III-  3) 


if  the  receiver  Y  has  local  oscillator  frequency  ta  ^  +  Aca^  anc*  Phase  Hence, 


-iAca  t  v 

X( CO )  V >••(-■)  x(w  -  coo)  y*(u>  -  coo>  e  0ellW'<P' 


co  >  0 


( III  —  4 ) 


In  order  to  obtain  the  eomplex  visibility  function  in  a  form  referenced  to  coordinates  fixed  in  the 
sky,  it  is  neeessary  to  take  out  the  rapid  phase  changes  eaused  by  souree  motion 


X(co)  Y*(co)  e+1‘<'Y  =  X(co )  Y*(co)  e  ref  ^  e1^'  'Y 


X(a  )  Y*(co)  e 


lea r  r  <ica/e)0  y  (ica/e)0  y 
ref  '  x'x  y  y 

e  e  J  J 


( III  —  5 ) 


i© 


The  visibility  function  referenced  to  the  sky,  A  e  ,  can  be  decomposed  into  various  frequency- 
dependent  and  independent  terms 


a  i®  v/  \  V/  w-  luTref 
A  e  =  X(ta)  Y  (a’)*  e 


id  t  c  i(ca-ca  )t  r  -lAta  t  v 

/  v  .  /  v  o  ref  o  ref  o  i(  0 - cp )  MTT  /  \ 

-  x(ca  —  d  )  y*(ca  —  d  )  e  e  e  e  v  (lil-b) 

oJ  o 


where 


ref 

d  — —  -  the  fringe  rotation  rate, 

{d  —  caQ)  t  j,  =  the  frequency  dependent  phase, 

Aa'o  =  the  local  oseillator  frequency  difference, 

(0  —  <p)  =  the  instrumental  phase. 

Equation  (III -6)  shows  that  it  is  possible  to  rotate  the  fringes  by  offsetting  the  local  oscillators, 
r  j,  should  also  inelude  constant  delays  or  instrumental  delays  in  the  system,  sueh  as  those  in 
cables  and  IF  amplifiers.  The  signals  are  filtered  so  that 

x1  ( co  —  d  )  =  x(ta  —  d  )  B  (ca  —  d  )  (III - 7 ) 

0  0X0 

y’(ca  -  dq)  =  y(D  -  caQ)  By(^  “  u Q)  ( III -8) 
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where  and  are  the  Fourier  transforms  of  the  bandpass  functions  for  the  two  receivers, 
so  that 


i© 

A  e  =  x’  (co  —  co  )  y*(co  —  co  ) 
'  o  J  o 


ico  t  r  1(0.1 -co  )t  _  -iAco  t  . 

o  ref  o  ref  o  i(G-<p) 

e _ e _ c _ e 

B  { co  —  co  )  B  *  ( co  —  co  ) 
x  o  y'  o 


(1H-9) 


This  may  be  computed  directly  from  a  Fourier  analysis  of  the  signals,  provided  the  time  varia¬ 
tion  of  the  phase  is  small  over  the  period  of  Fourier  analysis  T.  The  cross -spectral  function 
may  also  be  estimated  by  crosscorrelation  in  time  and  subsequent  transformation  since' 


(A  e1®}  ,  y 


+a 

a 


x'(t)  y'(t  -  Tref  -  T] 


ico  r  r  -iAco  t 

o  ref  e  o  ^i{Q-(p) 

B~1P 

x  y 


r  T 


- i( co-co  )t 

X  e  dT  w(t) 


(III -1  0) 


under  the  condition  that  T  »  a.  The  visibility  function  has  been  averaged  over  frequency  with 
an  effective  filter  response 


r: 


cos  Acot  w(t)  dT 


(III  - 1 1 ) 


where  w(t)  is  the  transform  weighting  function. 

The  time  crosscorrelation  technique  is  particularly  easy  to  implement  if  the  signals  are 
first  infinitely  clipped  so  that  when 


y‘,  x'(t)  >  0  yM,  x"(t)  -  1 
<0  yM,  x"(t)  =  -1 


(111-12) 


in  which  case 


<x’(t)  y’(t  -  T)>x  =  i<x2’(t)>T  <y2,(t)>Tl1  2 


x  sin  j  <xM(t)  yn(t  -  t)> 


T 


(111-13)’ 


from  the  theory  of  Van  Vleck.  The  time  crosscorrelation  can  only  be  performed  for  time  in¬ 
tervals  over  which 


ico  t  it)  -iAco  t 
o  ref  o 

e  e 


(III-14) 


However,  the  time  crosscorrelations  can  be  rotated  and  averaged  before  the  Fourier  transform 
is  taken.  In  practice,  the  local  oscillators  are  set  to  take  out  the  linear  portion  of  the  fringe 
rotation  and  the  rest  of  the  rotation  is  taken  out  after  crosscorrelation,  i.e.. 


x'(t)  y'(t  -  Tref  -  r)  e 


ICO  T  ~  -iAco  t 

ref  o 

c 


V 


x’(t)  y 


'(*  Trefr  T))t- 


lcor  it’)  -ico  t 
v  ref  o 

x  e  e 


(III  -1 5) 


t  x"  and  y"  are  normalized  to  unit  energy  for  time  T. 
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where  t?  is  a  discrete  time  midway  through  each  crosscorrelation.  An  additional  requirement 
for  this  scheme  to  work  is  that  the  change  of  from  one  period  to  the  next  is  small  enough 

that 

Au  Arref  «  7 r  .  (111-16) 


IV.  CALIBRATION  OF  AN  INTERFEROMETER 


The  normal  convention  in  interferometry  is  to  state  visibility  functions  normalized  so  that 
an  unresolved  point  source  has  a  fringe  amplitude  of  unity.  When  this  is  the  case,  the  brightness 
temperature  distribution  can  be  obtained  from  the  normalized  distribution  (the  transform  of 
normalized  fringe  amplitude  and  the  source  flux) 


Tat(u,  0,0) 
N  *  x'  y' 


-If 


A  (a: 
n 


(ia?/c)0  y  (Wc)9  y 

vx>  yy) e  e  y  y  dyxdry 


(IV-1) 


T„(w,  0,0)  = 
B  x  y 


F(co)  TNU',ex,  ey) 
/  TNdf2 


(IV-2) 


There  are  many  ways  of  obtaining  the  normalized  visibility  function,  about  the  simplest  of  which 
is  to  make  measurements  with  zero  spacing,  since 


A  ( a;,  y  , y  )  = 
n  '  x  '  y 


A(w,yx,yy) 
A  (u,  0,  0) 


( I V  —  3 ) 


but  it  is  seldom  possible  to  do  this  when  the  antennas  are  fixed.  Thus,  a  more  complicated 
method  of  calibration  has  to  be  employed  which  uses  the  measurements  made  with  the  individual 
antennas.  For  a  single  antenna 


ta<w>  = 


[x'(Cl)  —  CO  )  x 1  *  ( CJ  —  CO  ) 

1  o  o 


ref  ref 


X’ 

ref  ref 


(TV-4) 


where  T ^  is  the  system  temperature,  and  is  the  Fourier  component  of  the  video  output  on 

cold  sky.  From  Eq.(III-13), 


ta<“>  = 


Sn<-)TS-SnBp<-)TS 

nBP 


(IV —5) 


where  are  the  Fourier  transforms  of  the  clipped  autocorrelation  function 


s„i"'= i  “ini< 


xM(t)  xM(t  —  t)>t  e  1C^T  dr  w(r) 


and 


Tq  X  total  noise  power  on  source 

- 1  _  o _ 

S  total  noise  power  off  source 


When  the  projected  baseline  is  zero,  the  signals  are  perfectly  correlated  so  that 

JTAxlu)  Sx<")  +  J\n» 


x(  GO  ~~  00  ) 
O 


Jco) 


(IV-6) 


(IV-7) 
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( I V  —  8 ) 


y(co  —  co  )  =  c 
J  o'  y 


/ 


T.  (w)  S  (w)  +  /  T 

y  y 


ny(w)] 


where  c^,  are  constants,  and  s(co)  and  n(co)  are  normalized  so  that  they  each  have  unit  magni¬ 
tude.  From  Eqs.  (Ill —7 )  and  (III-8), 


so  that 


A(co,  0,  0)  =  c  c 

x  y 


/T A  (")TA 
NX  V 


N/(xrefx?ef>  ^ref^ref’  /TA  (w)  TA  (w) 
nx  y 


,/t  s  ts  yi 

N  x  y 


B  I2  |  B  |2 

xl  yl 


(IV  -9) 


ico  t  -  i(co-co  )t  f  -iAco  t  x 

,  ,  o  ref  o'  ref  o  i(0-<p) 

x’yT*  e  e  e  e 


A  ( co,  y  ,y  )  = 
n  x  ry 


/t„  t,  /b  2b  2 

/  S  S  v  x  y 
nx  y 


Sn  (W) 

*y 


B  B'1  /  T  A  tT  /Tx  rx  -  J  (y  -y  •'  -) 
x  y  /  A  An  ref  ref  JrefJref 
17  v  x  y 

/Ts  Ts 

N  x  y 


B  B  *  B  B* 
x  x  y  y 


BP 


‘BP 


(w)  (to)  Jta  (co)  ta  (co)  bvb* 


A 


A 


X  y 


(IV-10) 


(IV-11) 


for  the  clipped  system,  where 


Sn  <“> 
xy 


J  Sin  | 


x"(t)  y"(t-rref  ■ 

-i(co-co  )t 

x  e  °  dr  w(t) 


ICO  T  r  -lAcO  t  ~  x 
.  o  ref  o  i(0-<p) 

r  e  e  e 


(IV-12) 


V.  NOISE  ANALYSIS 

The  signals  can  be  represented  in  the  frequency  domain  by 


x(co)  = 

/ T.  (w)  S  (w)  + 

P7nx(aj) 

(V-l) 

J  X  N 

J  X 

y(«*>)  = 

J  Ta  («)  Sy(w)  + 

/rpy»> 

(V-2) 

where  and  are  the  antenna  and  system  temperatures,  respectively,  and  |S(co)|^  and 

i  T7 

|n(co)|  are  unity.  If  the  signal  and  noise  are  both  Gaussian,  they  have  independent  Gaussian 
components  with  probability  distribution 

,  r  2 

p(z)  =  -  e"z  (V-3) 

\flr 

where  z  =  Re  S(co),  Im  S(co),  Re  n(co),  Im  n(co)  so  that 
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(V-4) 


x(w)  x*(w)  =  T.  (w)  {[He  S  (to))2  +  [1m  S  (to))2} 

-rt.  X  X 


+  T„  {[Re  n  (to)]2  +  [Im  n  (to)]2} 

o  X  X 

X 


=  Ta  (to)  +  Ts 


while 


[ x( to )  x*(to)l2  -  (x(to)  x*(u>)]  =  T2  (to)  {[Re  S  (to)]4  +  [lm  S  (to)]' 


+  2  [Re  Sx(to)]2  [lm  Sx(to)]2}  4  T2  [Re  Sx(to)|4 


+  [lm  n^co)]  — 


ta  (w)  +  Ts 


4  2  [Re  Sx(o,')l2  [lm  n^co)]2 


4  4Ta  Ts  {[Re  Sx(to)]2  [Re  nx(to)|2  4.  .  .  } 
X  X 


=  Ta  (to)  4  T, 
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where  the  bar  denotes  the  statistical  or  ensemble  average.  If  an  average  is  performed  over  fre¬ 
quency,  then 


Var  <x(w)  x*(w)>Aw 
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+  T, 


A  CO 
Zi r 


At 


( V-6) 


since  there  are  (Aco/27t)  At  independent  frequency  Fourier  components  in  the  time  interval  At. 
The  above  result  is  the  well-known  noise  in  a  total-power  radiometer.  Now,  consider  the  cross 
spectrum  multiplied  by  a  factor  of  two  so  that  the  results  can  be  compared  with  a  total-power 
radiometer 


2x(w)  y*(«)  =  2  [TA  (w)  Ta  (w)  S  (w)  S*(o>)  +  2  / Tg  Tg  n  (w)  n*(w) 
v  x  y  ^  V  ‘  x  y  ^ 


+  2 


S  x'“"  y 

y  J 
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Sv(w)  n*( a; )  +  2  /  T^  (co)  TQ  nv(co)  S*(to) 


S  x 
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( V-7) 


The  first  term  is  the  signal  and  has  a  statistical  average  equal  to  twice  the  visibility  function, 
while  the  uncorrelated  portion  of  the  first  term  and  the  other  terms  make  up  the  noise.  The 
first  term  can  be  considered  as  the  vector 
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/ 


S  -  2  Re  / Ta  (a:)  Ta  («)  S  («)  S*(W)  i 
x  y  J  J 


+  2  Im  lTiX  ( Cc' )  Ta  (w)  S  (w)  S*(«)  i^ 
v  x  y  y 


In  the  case  of  a  weak  signal  (TA  «  T„)  perfectly  correlated  such  as  that  from  a  point  source, 
the  noise  vector  is  just 


N  -  2  He  /  Tc  Tc  n  (w)  n*(w)  i„ 
n/  Sx  Sy  X  y  Re 


+  2  Im  /  1  c  T~  n  (w)  n*(a>)  u 
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v  x  y 


( V-9) 


The  real  and  imaginary  components  are  Gaussian  random  variable  with  zero  mean  and 
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The  real  and  imaginary  components  are  independent  since 
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using  the  property  that 


X.X-X.X,  =  XjX-X.X.  +  X.X.X-X,  +  x.x.x.x, 

1234  1234  1324  1423 


(V-12) 


for  a  Gaussian  process.  The  probability  distribution  of  the  noise  vector  is 


,  ..  1  -(a2+b2),/2o2 

p(  a,  b)  =  - 2  e 

2na 


(V  -1  3) 


where 
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so  that 
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p(|Nl,=  1^1  e-|S|2/2»2  |5|>0 

(7 

P(0)  =  ^  0<  0  <  2tt  . 

The  length  of  the  noise  vector  is  a  Rayleigh  distribution  with  mean  \I  (ir/2)  a  and  (  |n|  2)  = 
shown  in  Fig.  1 . 


( V-l 7) 

( V  -1 8) 
o  2 

2cr  ,  as 


Fig.  1.  Probabil i ty  distribution 
of  amplitude  of  noise  vector. 


The  rms  deviation  of  the  amplitude  A  and  phase  (p  of  the  visibility  function  estimate  can 
be  simply  written  only  when  the  signal  vector  is  much  stronger  than  the  noise  vector,  in  which 
case 

|s  +  N|  ~\s\  ( V -1 9) 

|s  +  N|2  =  |s|2  +  |n|2  ( V  -20) 


and 


AA  =  (|N|2)l/2=  4T4  = 
rms  1  1 
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A  =  2  /TA  (c o)  T^  (c o)  °K  for  point  source 
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(V-22) 


and 


A  (p 


rms 


|N|  2  sin2  0 

X2 


42 


U 

A 


[AuT 
\l  2rr 


( V-23) 


At 


Expressions  for  intermediate  signal  levels  have  to  be  evaluated  by  integration.  If  the  cross¬ 
correlation  technique  is  applied  to  a  single  antenna,  the  signal  vector  reduces  to  /  T  A  T  A  since 

v  Ax  Ax 

3  dB  is  lost  in  splitting  the  power  between  receivers,  while  the  rms  noise  is  a  =  (Aco/27t)  At 

since  only  one  component  of  noise  affects  the  signal  whose  phase  is  fixed.  The  crosscorrelation 
receiver  thus  has  \  2  more  noise  than  the  total-power  receiver,  but  yields  the  additional 
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information  of  receiver  baseline  subtraction.  When  an  interferometer  is  used  to  measure  the 
flux  of  a  point  source  whose  position  in  the  fringe  pattern  is  unknown,  the  signal-to-noise  ratio 
is  the  same  as  that  of  a  total-power  radiometer  of  one  antenna  if  both  interferometric  antennas 
and  receiver  systems  have  the  same  signal-to-noise  ratios.  In  the  preceding  analysis,  only  the 
case  of  a  point  source  and  a  high  system-to-signal  temperature  ratio  was  considered.  In  the 
other  extreme,  when  the  normalized  fringe  amplitude  is  much  smaller  than  unity,  the  magnitude 
of  the  signal  vector  is 


2  IT  (W)TA  («j  I A  | 

NX  y 


and  the  rms  noise 
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At 


(V-24) 


( V  -25) 


since  the  uncorrelated  signals  now  effectively  add  to  the  system  temperatures.  This  result  can 
be  seen  by  evaluating  the  cross  terms  of  Eq.  (V-7)  which  were  previously  neglected. 


VI.  SOURCE  MODELS 

Certain  source  distributions  have  very  simple  complex  visibility  functions.  For  example, 
a  point  source  has  a  normalized  visibility  magnitude  of  unity  and  a  phase 

r-^  {( sin  6^  cos  6^)  A6^  +  [cos6^  cos6s  sin(L^  -  LR)]  ARAS 

-  [cos6r  sin  6^  cos(Lg  -  Lg)]  A6g}  (VI-1) 

where  ARA<.  and  A6^  are  the  position  offsets  of  the  source  from  the  reference  position.  Even 
if  only  a  small  coverage  of  projected  baseline  is  observable,  it  is  generally  possible  to  use  a 
least-squares  fitting  technique  to  find  the  source  position.  If  the  signal-to-noise  is  good,  no 
ambiguities  are  likely  to  arise  in  the  fitting  process,  especially  if  the  phase  is  tracked  continu¬ 
ously  so  that  only  a  small  range  of  multiples  of  Zn  need  be  tried  in  the  fitting  procedure.  In 
general,  the  fitting  process  is  better  than  taking  the  two-dimensional  transform  only  over  the 
range  of  projected  baselines  covered,  owing  to  the  restriction  to  a  particular  source  model. 

Other  simple  visibility  functions  are  those  of  a  circularly  symmetric  source  whose  phase 
is  zero  when  refined  to  the  source  centroid,  and  a  uniform  disk  source  of  radius  R  whose  fringe 
amplitude  is 

2tt) 

E 1 

S 

where  is  the  first-order  Bessel  function,  and  S  is  the  fringe  spacing  or  the  reciprocal  of  the 
projected  baseline  length  in  wavelengths. 

VII.  INDEPENDENT  TIME  STANDARDS  INTERFEROMETER 

Recently,  several  interferometers  have  been  operated  with  independent  time  and  frequency 
standards  at  the  two  antennas,  thereby  eliminating  the  real-time  link  between  the  two  antennas. 
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Coherent  integration  of  the  visibility  function  over  a  time  T  is  possible  if  the  spectral  purity 
and  stability  of  the  two  frequency  standards  are  such  that 

!  <eiA<p(t))Tj  ^  l  (V1I-1) 

where  Ac  (t)  is  the  phase  difference  between  the  two  local  oscillators  after  removal  of  any  con¬ 
stant  frequency  difference  by  a  frequency  offset  search.  For  a  rubidium-controlled  crystal  clock 
(a  crystal  clock  which  is  phase  locked  to  a  rubidium  line),  coherent  integration  times,  for  a 
system  at  1666  MHz,  of  a  few  minutes  are  possible.  For  a  resolution  bandwidth  of  Af,  the  max¬ 
imum  allowable  time  error  At  must  be  such  that 

AfAr  «  1  .  (VII -2) 

I'ime  synchronization  close  enough  for  an  initial  signal  detection  can,  in  practice,  be  achieved 
using  Horan  time  transmissions.  A  fringe  frequency  and  time  error  search  can  be  performed 
by  computing  the  visibility  functions  for  short-time  intervals  with  sufficient  frequency  resolution 
to  account  for  any  delay  errors  (or  equivalently  taking  enough  points  in  the  crosscorrelation 
function  to  ensure  reaching  the  real  zero  point),  and  then  coherently  averaging  the  visibility 
funetion  to  maximize  the  amplitude  and  minimize  the  phase  shift  across  the  band. 

Stated  mathematically,  the  frequency  offset  and  delay  are  determined  by  maximizing 


l  Q.Aw)  . 

V  V  \  /x  t,V  lb  t  IT  Cl' 

L,  \t,(w)  e  e  e 

t'  U? 


( VII  -  3) 


iO  f(w) 

where  A  ^ ,  ( u: )  e  is  the  complex  fringe  amplitude  determined  over  a  short  period  centered 

at  t’.  t  and  ]  are  the  delay  and  frequency  search  parameters,  which  may  be  quantized  provided 

|  <e1<A,'t/2»)>T|  w  t  (VII -4) 

and  AtAu.  «  7 r.  The  delay  search  is  more  conveniently  performed  in  the  time  domain  by  search¬ 
ing  the  delay  shift  of  the  crosscorrelation  function. 

For  a  continuum  source,  the  cross  spectral  function  of  the  video  signals  is 


...  27Tll<  t  1WT.  1C x3T  -IU)T  . 

.  it)  K  intf  e  st 

A  e  e  e  c  e 


(VII-5) 


for  the  single -sideband  system  and 


io;t.  ,r  leer  -icjt  , 
.  ,  o  n*  *\  intf  c  si 

A  cos  0  +  27rF^t)  e  e  e 


(Vll-6) 


for  a  double- sideband  system,  where  \  and  O  are  fringe  amplitude  and  phase,  F^  is  the  fringe' 
rate,  t.  ^  is  the  geometric  interferometer  delay,  and  t^  is  a  timing  error,  t^  is  the  tape  time 
shift  necessary  to  take  out  T^ntf  The  double- sideband  system  has  the  disadvantage  of  requiring 
a  least-squares  fitting  procedure  owing  to  the  lack  of  a  quadrature  component,  but  has  the  ad¬ 
vantage  of  timing  errors  not  affecting  the  fringe  phase.  The  Fourier  transform  of  ( V 1 1  —  5 )  is 

R  (t)  =  <x(t)  y(t  -  r)>  2A  cos ( O  +  2?rFRt)  +  T.ntf  +  rp  -  Tgt) 

-  2A  sin(Q  +  27rFRt)  g2(r  +  Tjntf  +  Tp  -  rgt)  ( \  I1  -7) 
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where 


and 


■  I 


;  ,  sin  a'  t 

max  _  (ioj  max 

cosier  -=—  =  - - - 

2ir  2t tt 


.  *  f^rnax  .  die 

g2(r)=j  Sln^  -27 


1  —  COS  Cl  r 

max 

2  7T  T 


( Vll-8 ) 


( \  11-9) 


while  the  transform  of  (VII -6)  is 

Kxy(T)  -  A  eos(0  +  27rFRt)  g^r  +  +  tq-  t^)  .  (VI1-10) 

The  fringe  rate  FR  includes  both  the  source  motion  and  the  oscillator  offsets  and  is  sufficiently 
slow  that  immediate  erosseorrelation  may  be  performed.  Additional  rotation  and  delay  are  intro¬ 
duced  in  the  fringe  rate  and  delay  search  as  in  Eq.(VIl-3),  as  often  as  is  required  to  ensure  no 
significant  reduction  in  fringe  amplitude. 

The  interval  T  over  which 


<eiA</,(t)>  |*1 


(V1I-11) 


is  the  maximum  coherent  integration  period  possible  without  significant  reduction  in  fringe  ampli¬ 
tude.  After  the  limit  of  coherent  integration  has  been  reached,  it  is  still  possible  to  improve 
the  knowledge  of  the  fringe  amplitude  with  incoherent  averaging.  The  incoherent  average 
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yields  no  phase  information  and  cannot  be  used  for  position  measurements  or  to  obtain  a  bright¬ 
ness  distribution,  but  the  magnitude  will  still  indicate  the  effective  source  size.  The  signal-to- 
noise  ratio  that  can  be  obtained  with  incoherent  averaging  approaches  that  of  the  Hanbury  Brown 
and  Twiss  interferometer  in  the  limit  when  the  bandwidth  of  the  coherent  integration  B\V  ap¬ 
proaches  the  bandwidth  of  signals  before  erosseorrelation  —  in  other  words,  when  the  coherent 
integration  time  reduces  to  the  reciprocal  of  the  pre-erosseorrelation  bandwidth. 

For  a  spectral  line  source,  the  noise  analysis  is: 
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from  Fqs.(V-20)  and  (V-21).  For  small  signals. 
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If  BW  is  defined  as  the  reciprocal  of  At,  and  t  as  the  total  integration  time,  then  FCq.  (Vll-14) 
becomes 
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and  the  signal-to-noise  ratio  is 
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When  the  source  being  observed  has  spectral  features,  it  is  possible  to  coherently  average 
the  data  using  one  of  the  strong  features  as  a  phase  calibration.  This  method  can  effectively 
make  the  interferometer  phase  stable  when  the  signal-to-noise  ratio  is  good  on  the  calibration 
feature.  However,  the  eomplex  visibility  functions  are  now  no  longer  referred  to  the  sky,  and 
a  map  relative  to  the  calibration  feature  will  be  obtained.  For  this  averaging  process. 
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where  ^  ^  is  the  phase  of  the  calibration  feature  for  the  n^  coherent  integration  period.  This 

process  has  no  effect  on  the  interferometer  noise  provided  the  calibration  feature  has  negligible 
noise.  Another  possible  technique  for  removing  differential  phase  noise  due  to  drifting  standards 
is  to  simultaneously  observe  a  strong  calibration  souree  with  another  set  of  antennas  or  off-axis 
feed  system.  If  the  calibration  source  is  strong  enough  to  produce  a  phase  with  little  noise  in 
the  coherent  integration  period,  these  phases  ean  then  be  used  to  increase  coherent  integration 
time  on  the  souree  being  examined.  Of  eourse,  the  different  position  of  the  calibrator  would 
have  to  be  taken  into  account  in  obtaining  the  phase  correction  term.  This  technique  could  also 
be  used  to  remove  the  atmospherie  distortion  which  limits  optical  astronomy.  For  example,  if 
a  calibration  star  is  sufficiently  close  to  the  objeet  of  interest,  say,  a  planet,  the  fringe  pattern 
will  be  jittering  in  the  same  fashion  on  the  planet  as  it  does  on  the  star;  consequently,  it  should 
be  possible  by  electronic  image  superposition  to  integrate  the  visibility  functions  coherently. 

The  phase  ehange  due  to  the  atmosphere  to  the  first  order  for  a  plane  parallel  atmosphere  is 


(jj 

*e 


rL 

rL'  1 

\  (n(i  )  -  1)  d /  - 

\  (n(f)  -  1)  di'} 

Jo 

Jo  | 

(V1I-16) 


where  f  and  f  are  the  paths  through  the  atmosphere.  Typieally,  the  atmosphere  adds  about 
20  cm  to  the  path  and,  owing  to  turbulenee,  the  path  will  not  be  the  same  for  the  two  apertures 
of  the  interferometer.  This  differential  phase  fluctuation  whieh  is  given  by  the  above  equation 
is  not  serious  at  radio  frequencies  below  about  3  GHz  as  it  seldom  reaches  90°  peak,  even  for 
widely  separated  antennas.  However,  at  optieal  frequencies,  it  is  significant  for  paths  sep¬ 
arated  by  only  a  few  eentimeters.  An  approximate  angular  distance  over  which  the  differential 
phase  fluctuations  will  remain  highly  correlated  is  the  diameter  of  an  aperture  over  which  phase 
deviations  are  smaller  than  ir  divided  by  the  length  of  the  atmosphere  where  most  fluctuations 
oeeur.  In  the  optieal  range,  this  is  only  about  20  seeonds  of  are,  but  this  is  still  large  compared 
with  the  resolution;  consequently,  the  use  of  a  phase  calibrator  within  20  seconds  of  the  object 
of  interest  eould  be  used  to  take  out  atmospherie  fluctuations.  At  8  GHz,  a  phase  calibrator 
within  30  minutes  of  are  would  probably  suffiee. 
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VIII.  HAYSTACK-MILLSTONE  INTERFEROMETER 


The  120-foot  antenna  of  the  Haystack  Microwave  Research  Facility  and  the  84-foot  antenna 
of  the  Millstone  Radar  Facility,  both  operated  by  Lincoln  Laboratory,  were  used  as  an  inter¬ 
ferometer.  The  antennas  are  separated  by  approximately  2250  feet  along  a  line  approximately 
19°  East  of  North.  This  baseline  gives  a  minimum  fringe  spacing  of  54  seconds  of  arc  at  18  cm 
and  provides  a  good  range  of  projected  baseline  for  a  wide  range  of  declination.  The  convenient 
baseline  and  the  spectral  processing  equipment  at  Haystack  make  the  system  ideal  for  a  study 
of  OH  emission  regions  that  were  unresolved  with  a  single  antenna.  Figure  2  shows  the  projected 
baseline  coverage  for  the  OH  emission  sources. 


Fig.  2.  Fringe  amplitude  phase  coveroge 
for  some  OH  sources.  Limits  shown  are 
observoble  limits. 
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A  block  diagram  of  the  Millstone  receiver  front  end  is  shown  in  Fig.  3.  Circular  and  linear 
polarization  is  obtained  from  a  dual-mode  horn  whose  output  ports  give  vertical  and  horizontal 
polarization.  Combining  vertical  and  horizontal  signals  through  a  hybrid  complex  gives  left- 
and  right-circular  polarization  after  correct  adjustment  of  the  phase  lengths.  The  antenna  output 
is  amplified  by  a  tunnel  diode  amplifier  and  then  filtered  to  reject  the  image  band.  A  ferrite 
switch  is  included  for  calibration  measurements.  In  the  normal  interferometer  mode,  the  switch 
remains  switched  to  the  antenna  side.  A  noise  source  is  used  for  single-antenna  and  system- 
temperature  measurements. 

The  local  oscillator  signal  is  derived  by  phase  locking  an  oscillator  to  the  sum  or  difference 
of  a  harmonic  of  67  MHz  and  a  signal  whose  frequency  could  be  varied  approximately  28  MHz. 

For  the  OH  emission  measurements,  the  24  harmonic  was  selected.  The  output  of  the  oscillator 
was  filtered  to  attenuate  any  spurious  signals  that  tend  to  be  produced  by  the  synchronizer. 

Figure  4  is  a  block  diagram  of  the  local  oscillator  system. 

The  mixer  output  is  amplified  by  a  30-MHz  amplifier  with  a  10-MHz  bandwidth.  A  line  driver 
then  boosts  the  level  to  lOOmW.  The  Haystack  receiver  front  end  did  not  require  the  line  driver, 
owing  to  its  proximity  to  the  control  room  where  the  IF  outputs  are  combined.  Otherwise,  the 
Haystack  front  end  is  similar  to  that  at  Millstone. 

The  intersite  coupling  of  the  radiometers  involves  the  transmission  of  antenna-pointing 
commands  to  Millstone,  remote  control  of  the  radiometer,  transmission  of  the  intermediate 
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Fig.  3.  Radio  frequency  section  of  Millstone  receiver. 


Fig.  4.  Locol  oscillotor  system  shown  for  center  of  band  ot  1666MHz. 
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Fig.  5,  Line  serva  to  maintain  phase  stability. 

frequency  to  Haystack,  and  two-way  coupling  of  the  local  oscillator  reference  signals.  To  main¬ 
tain  good  phase  stability,  the  reference  signals  are  transmitted  along  a  servo- controlled  line 
that  nullifies  line  length  changes  caused  by  temperature  change  and  other  effects.  The  trans¬ 
mission  system  also  has  to  overcome  a  large  line  attenuation  of  55  dB  for  one-way  transmission 
at  67  MHz.  The  selection  of  a  lower  basic  frequency  would  have  reduced  attenuation,  but  would 
have  made  phase  locking  to  L-band  more  difficult.  A  block  diagram  of  the  line  servo  is  shown 
in  Fig.  5.  The  67-MHz  reference  signal  is  amplified  to  approximately  1  watt  and  transmitted  to 
the  line  through  a  hybrid  junction.  At  the  receiving  end  of  the  line,  a  portion  of  the  signal  is 
reflected.  The  reflected  signal  undergoes  phase  reversal  with  a  100-kHz  rate  as  the  diode  switch 
modulates  the  reflection  coefficient  from  to  — Very  little  of  the  100-MHz  modulation  is 
passed  into  the  67-MHz  amplifier  because  of  the  isolation  afforded  by  the  hybrid  tee.  At  the 
transmitting  end,  the  reflected  signal  is  mixed  with  the  67-MHz  reference  signal  and  amplified. 
The  output  of  the  100-kHz  amplifier  is  proportional  to  cos  q  cos  27rft,  where  f  =  100  kHz,  and  <p 
is  the  phase  of  the  67-MHz  signal  after  having  traveled  twice  the  line  length.  Multiplication  of 
this  signal  by  the  100-kHz  reference  signal  produces  the  necessary  error  signal  from  the  servo 
loop.  The  high  loop  gain  makes  it  possible  for  the  line  servo  to  maintain  a  constant  electrical 
line  length  to  within  a  small  fraction  of  an  inch. 

The  IF  signals  are  either  added  or  subtracted,  as  shown  in  Fig.  6.  The  combined  signal  is 
filtered  and  converted  to  video.  The  autocorrelation  function  of  the  clipped  and  sampled  video 
signal  is  taken  with  a  digital  correlator. 

The  difference  of  the  autocorrelation  functions  for  the  added  and  subtracted  signals  yields 
the  real  part  of  the  crosscorrelation  function  since 

[ x( t )  +  y(t)|  [x(t  -  t)  +  y(t  -  t )]  -  [ x( t )  -  y( t )]  [x(t  -  r)  -  y(t  -  t)] 

=  2x(t)  y(t  -  t)  +  2y(t)  x(t  -  t)  (VIII-1) 

and 
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Fig.  6.  Crosscarrelotian  scheme  by  addition  ond  subtraction  of  IF  signols. 


Fig.  7.  Hoystock-Mil Istane  interferometer. 
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2  \  (<x(t)  y(t-r)>T  +  <y(t)  x(t-r)>T]  e'la;T  dr 
=  4  Re  x(w)  y*(w) 

*iuJTref(t) 

-  4  Re  A  e  1  .  (VIII -2) 

The  real  and  imaginary  components  of  the  eomplex  visibility  function  can  be  found  by  a  least- 
squares  fitting.  This  system  yields  an  effective  system  temperature  of 

v  x  y 

if  the  signal  powers  are  added  with  equal  power.  However,  the  noise  level  is  twice  that  of  direct 
erossc  or  relation. 

Autocorrelation  functions  of  the  sum  and  difference  signals  are  transferred  alternately  to 
a  computer  every  100msec.  The  basic  correlation  period  is  87.5  msec.  The  system  is  blanked 
for  12.5  msec  while  the  data  are  being  transferred  to  the  U490  computer.  Signal  bandwidths  of 
4  MIIz,  1.2MIIzf  400  kHz,  120  kHz,  and  40  kHz  can  be  analyzed  by  selecting  appropriate  filter 
and  sampling  frequencies.  The  autocorrelation  function  is  a  1 6-bit  binary  word  for  each  of  100 
delays.  An  extra  bit  is  used  to  indicate  the  state  of  the  switching  reference  signal.  The  auto¬ 
correlation  functions,  together  with  continuum  data,  bandwidth  and  correlator  mode,  and  antenna 
command  azimuth  and  elevation,  are  transferred  to  magnetic  tape.  A  block  diagram  of  the  whole 
system  is  shown  in  Fig.  7. 

Single-antenna  measurements  were  made  to  measure  the  aperture  effieieneies.  The  antenna 
temperatures  of  Cassiopeia  A  were  160°  and  180°K,  which  indicated  efficiencies  of  40  and  25  per¬ 
cent  for  the  Millstone  and  Haystack  antennas,  respectively. 

After  the  individual  radiometers  were  examined  for  linearity,  bandpass,  and  lack  of  spurious 
signals,  the  system  was  checked  for  phase  coherence  and  phase  stability.  The  test  arrangement 
is  illustrated  in  Fig.  8.  A  signal  from  a  test  oscillator  is  connected  to  each  radiometer  individ¬ 
ually  and  then  simultaneously.  When  the  oscillator  is  swept  in  frequency,  the  bandpass  is  dis¬ 
played  when  only  one  radiometer  is  connected  to  the  oscillator.  When  both  radiometers  are 


Fig.  8.  Interferometer  test  system. 
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swept  simultaneously,  fringes  appear  because  of  the  phase  relationship  between  the  radiometer 
outputs.  The  fringe  spacing  in  Hertz  is  the  reciprocal  of  the  delay  in  seconds.  Measurements 
of  the  fringe  pattern  and  bandpass  functions  showed  that  the  phase  noise  in  the  system  produced 
less  than  2-percent  reduction  in  fringe  amplitudes.  The  phase  stability  was  better  than  50°  during 
24  hours.  Some  of  this  shift  may  have  been  due  to  the  measuring  technique,  as  observations  of 
continuum  radio  sources  indicated  better  stability. 

The  trombone  section  of  the  line  servo  was  observed  to  move  approximately  5  feet  during 
the  sunrise  and  sunset  period  when  the  temperature  changed  about  50°  F. 


IX.  WIDE  EFFECTIVE  BANDWIDTH  INTERFEROMETER 

When  an  independent  standards  long-baseline  interferometer  is  used  to  measure  the  positions 
of  unresolved  continuum  point  sources,  the  phase  information  is  seldom  usable  since  it  is  not 
possible  to  calibrate  the  phase  at  time  intervals  sufficiently  frequent  to  ensure  its  constancy. 
Consequently,  the  position  must  be  obtained  by  least-squares  fitting  the  fringe  rate  F^  to  obtain 
the  offsets  in  right  ascension  and  declination.  From  Eq.(II-l6), 


dwr. 


fr  = 


intf 


dt 


(  ,  T  }  ^ 

—  “dT  cos6b  cos6s  sin(Ls"LB} 


( I X  - 1 ) 


and,  hence. 


AF 


R 


WD  S  r-  r  A  A  •  r  .  ,  . 

—  — tt  coso,,  [Ao0  sino(^  siniLc 
c  dt  B  1  b  S  S 

+  ARAS  cos  cos(  Lg  —  LR)] 


L„) 


( I X  —  2 ) 


If,  however,  the  delay  can  also  be  precisely  measured,  then  least-squares  fitting  to  the  delay 
difference  At^^.  yields  more  information  on  the  source  offsets.  From  Eq.  (II -16), 


At 


intf 


i-  {(sin <5^  cos  6^)  A6g  +  [cos  6^  cos  6^  sin ( ,Lg  —  L^)j  ARAg 
—  [cos  6  sin  6S  cos  (  Lg  —  L^)]  A6g} 


( I X-  3) 


Like  the  fringe-rate  fitting,  the  delay  does  not  have  the  ambiguity  difficulties  that  fringe-phase 
fitting  presents.  In  fact,  one  measurement  of  fringe  and  delay  gives  the  source  position  offset. 
Fringe-rate  fitting  yields  an  offset  resolution  of 


X. 

D 


AFr 


27r(dLs/dt)J 


(IX— 1) 


or,  if  the  fringe  rate  can  be  measured  to  1  cycle  per  hour,  the  offsets  can  be  measured  to  ap¬ 
proximately  24  times  the  fringe  spacing.  Delay  fitting  yields  an  offset  resolution  of 


where  A. and  a; o/ 2 7r  are  some  "center"  wavelength  and  frequency.  For  example,  measurements 
of  delay  to  1  nsec  would  yield  an  offset  resolution  of  twice  the  fringe  spacing  for  a  2-Gllz  center- 
frequency  system. 

In  practice,  measurement  of  delay  between  two  very  long  baseline  interferometer  stations 
would  not  require  recording  a  very  large  bandwidth,  as  the  signal  can  be  sampled  at  various 
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small-bandwidth  windows  within  a  large  bandwidth.  Extrapolation  of  the  fringe  phase  from  one 
window  to  the  next  could  provide  an  increasingly  precise  delay  measurement. 

Measurement  of  delay  and  fringe  rate  with  a  very  long  baseline  interferometer  would  provide 
an  extremely  precise  technique  for  measurement  of  the  positions  of  radio  sources;  then,  re¬ 
versing  the  previous  analysis  would  provide  a  precise  technique  for  measurement  of  positions 
on  the  earth's  surface  as  well  as  irregularities  in  the  earth's  rotation. 

X.  SUMMARY  AND  CONCLUSIONS 

In  Sec.  II,  we  discussed  the  principles  of  interferometry,  or  how  the  angular  distribution  of 
electromagnetic  radiation  is  related  to  the  electric  field  at  certain  "sample"  points.  While  an 
antenna  samples  and  adds  the  electric  field  at  many  points  by  virtue  of  its  structure  or  geometry, 
a  two-element  interferometer  samples  the  single  antenna  "filtered"  electric  field  at  only  two 
points  instantaneously,  but  can  build  up  many  sample  points  as  the  baseline  changes  in  time. 

Since  all  the  sample  points  are  not  merely  added  but  crosscorrelated,  the  resultant  antenna  can 
be  imagined  to  be  electrically  steered  by  the  phase  term  in  the  two-dimensional  transform. 

The  noise  analysis  of  Sec.  V  shows  the  relation  of  the  noise  in  the  fringe  amplitude  and  phase 
to  the  noise  levels  of  the  individual  systems.  It  is  interesting  to  note  that,  while  an  interferometer 
can  map  an  area  within  the  single-antenna  beam  with  a  signal-to-noise  ratio  of  a  single  antenna, 
a  super  large  parabolic  antenna  (diameter  approximately  that  of  the  maximum  interferometer 
baseline)  used  to  map  the  same  region  does  not  gain  in  signal-to-noise  ratio  by  the  ratio  of  its 
area  to  the  single-dish  area.  This  is  because  the  super  large  dish  has  to  be  mechanically 
scanned  over  the  region;  hence,  the  effective  integration  time  is  reduced  by  the  ratio  of  its  beam 
area  to  the  map  area.  Thus,  an  interferometer  uses  the  collecting  area  efficiently,  like  a 
multiple-feed  antenna  or  camera  with  photographic  emulsion. 

Although  this  report  mainly  discusses  spectral  line  interferometry  or  the  mapping  of  fre¬ 
quency  dependent  brightness  distributions,  the  section  on  wide  effective  bandwidth  interferometry 
(Sec.  IX)  can  be  alternately  viewed  as  using  the  wide  bandwidth  to  increase  the  number  of  sample 
points  in  the  complex  visibility  plane  when  the  source  distribution  remains  constant  with 
frequency. 

Because  very  large  antenna  structures  are  not  mechanically  feasible,  interferometry  has 
become  a  very  important  technique  for  high  resolution  source  mapping  and  point  source  position 
measurements.  Very  long  baseline  interferometers  have  application  for  clock  synchronization, 
antenna  site  position  measurement,  and  numerous  geophysical  experiments  as  well  as  astronomi¬ 
cal  measurements. 
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